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COMPARISON OF ABSOLUTE CALCULATED AND MEASURED GAMMA AND NEUTRON 
DOSES IN TUNGSTEN - WATER-MODERATED CRITICAL ASSEMBLY 
by Paul G. Klann and Walter A. Paulson 
Lewis Research Center 

SUMMARY 

The mixed gamma and neutron dose was measured at 60 locations within the 
Tungsten - Water -Moderated Reactor critical assembly. Graphite wall thimble ioniza- 
tion chambers filled with carbon dioxide and polyethylene wall chambers filled with ethyl- 
ene were used. The chambers were absolutely calibrated in a bremsstrahlung beam 
against a secondary calibration standard and in a reactor against a water filled calorim- 
eter. These calibrations were used to obtain gamma and neutron response coefficients 
for the chambers which permitted partitioning of the measured mixed radiation dose into 
a gamma dose and a neutron dose. The measured gamma doses were compared with an 
ATHENA Monte Carlo calculation. Good agreement was found for the 16 locations com- 
pared. The average deviation was 10 percent. In addition, the measured neutron doses 
were compared with a "first collision" calculation of the dose. The calculated neutron 
doses were uniformly low with an average deviation from the measurement of 18 percent. 


INTRODUCTION 

One of the important aspects of the Tungsten - Water-Moderated Rocket Reactor fea- 
sibility program was the specification of the gamma dose in the core water moderator and 
the core structure. A Monte Carlo gamma heating digital computer program, ATHENA, 
was developed by the United Nuclear Corporation for the NASA Lewis Research Center to 
calculate the rate of gamma energy deposition as a function of position in the core. 

The accuracy of this program and the cross section data were checked by comparing 
gamma dose measurements conducted by the General Atomic Division of General Dynam- 
ics at the Tungsten - Water-Moderated Reactor (TWMR) critical facility with values cal- 
culated using ATHENA. The comparison of the calculated and measured gamma dose is 
given in this paper. 



In addition, measurements were conducted by General Atomic in which the total dose 
from combined gamma and neutron radiation was determined at a number of locations 
within the critical assembly. The neutron dose results are compared with a refined first 
collision dose calculation in this paper. 


GAMMA DOSE CALCULATIONS 
Description of TWMR 


Sixteen graphite wall, carbon dioxide filled thimble ionization chambers were used to 
measure the absolute gamma dose at 60 locations within the control tubes and the fuel 



Figure 1. - Cross section of TWMR critical assembly. 
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elements of the Tungsten - Water-Moderated Reactor (TWMR) critical assembly (refs. 1 
and 2). The TWMR critical assembly consisted of a hexagonal array of 121, 2. 562-inch- 
(6. 507 cm) diameter, 42. 5-inch- (108 cm) long fuel elements spaced on a 3-inch 
(7.62 cm) triangular pitch in the water moderator. The core was reflected by beryllium 
backed with water on all sides with exception of the top 2.0-inch (5. 08 cm) water reflec- 
tor. A cross section of the core is shown on figure 1. The bottom beryllium reflector 
which was built up of 2g- by 2^ by 4 inches (7.30 by 7. 30 by 10. 16 cm) beryllium blocks 
placed in the space between the core support plate and the grid plate is not shown on this 
figure . 

Each fuel element was composed of a series of concentric uranium 235, tungsten, 
aluminum, and uranium 238 rings within a hollow water-tight aluminum support tube. 
Figure 2 shows the cross section of the fuel element. Control tubes containing dilute so- 
lutions of cadmium nitrate were located in the triflutes between the fuel elements. 


42.5 in. 
(108 cm) 


70 in. 
(175 cm) 


Subassembly 
center post^ 



Fuel spacer tube - 

Tungsten tube 
(0.005 wall) — 


Subassembly 
lower spacer tube- 



'Turn-tite' 

closure 

-Subassembly upper 
spacer tube 

Subassembly 
end plate 

/-Pressure 
/ tube 

---Zirconium 

stud 

" — Zirconium 
center divider 


■Spacer ring 

Subassembly 
end plate 


Zirconium 
lower plate 


Pressure tube 
lower end cap 


2.562 in. 
(6.51 cm) 


Figure 2. - TWMR fuel element. 
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Computer Program Description 


The digital computer program, ATHENA, was used to calculate the gamma energy 
deposition rates at 16 of the graphite thimble ionization chamber locations. This com- 
puter program (ref. 3) employs Monte Carlo techniques to compute gamma-energy depo- 
sition rates as well as neutron and gamma fluxes in selected energy groups in complex, 
three-dimensional geometries. The option of representing a reactor by a symmetric 30° 
sector of the core with reflective boundaries at 0° and 30° makes this program especially 
efficient for calculations involving hexagonal-lattice reactors such as the TWMR critical 
assembly. 

In ATHENA, the source neutrons and photons are selected based on detailed input- 
specified reactor geometry, region and importance sampling parameters, power pattern, 
and operating history. These are followed through the various interactions until the oc- 
currence of capture, death by energy or region importance, or a degradation below an 
energy cutoff. 


Calculationa! Model 

The TWMR critical assembly was a heterogeneous core. Although the storage limi- 
tations of the digital computer did not permit detailed representation of each fuel cylinder 
in a stage, the essential heterogeneous nature of the core was retained in the calcula- 
tional model. The main features of the core representation are as follows: 

(1) The critical assembly was represented as a 30° sector of the hexagonal core with 
reflective boundaries at the 0° and 30° planes. 

(2) Each fuel assembly was represented as a single right cylinder having the homog- 
enized fuel element composition. 

(3) The aluminum pressure tubes, water moderator, and bottom and side reflectors 
were represented explicitly. 

(4) Each poison control tube was represented as a homogenized cylinder of alumi- 
num plus cadmium control solution. 

(5) Sixteen thimble ionization chambers were represented as right cylinders. Pre- 
liminary calculations showed that the small active volume of the chambers (about 1.8 cc) 
did not accumulate a good statistical sample of track lengths during the Monte Carlo cal- 
culation. To improve the statistics, the representation of the active volume was in- 
creased to 8.9 cubic centimeters by increasing the length. The experimentally mea- 
sured axial distribution of absorbed dose is a smooth function of axial location; hence, 
the increased length of the active section of the chamber should not introduce a signifi- 
cant error in the calculated dose. 
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TABLE I. - COMPOSITION USED IN MONTE CARLO 


CALCULATIONS 



Isotope 

Atom density 



atom/barn-cm 

atom/c c 

Fuel element 

Uranium 235 
Uranium 238 
Tungsten 182 
Tungsten 183 
Tungsten 184 
Tungsten 186 
Aluminum 
Zirconium 
Nickel 

7.5992X10' 4 

3.1228X10' 3 

5.1651X10' 4 

2.8173X10' 4 

5.9869X10" 4 

5.5839X10' 4 

2.8468X10' 2 

1.2695XKT 4 

8.8814X10' 5 

7.5992X10 20 
3.1228X10 21 
5.1651X10 20 
2.8173X10 20 
5.9869X10 20 
5. 5839X10 20 
2.8468X10 22 
1.2695X10 20 
8. 8814X10 19 

Poison tube 

Hydrogen 

Oxygen 

Aluminum 

Cadmium 

5.275 xlO' 2 
2.6375XKT 2 
1.26 xlO' 2 
1.1282X10" 4 

5.275 XlO 22 
2.6375X10 22 
1.26 XlO 22 
1. 12 82x1 0 20 

Dosimeters 

Carbon 

6.1889X10' 2 

6.1889X10 22 

Polyethylene 
and poly- 
styrene 

Carbon 

Hydrogen 

3.99 xlO' 2 
7.98 xlO" 2 

3.99 XlO 22 
7.98 XlO 22 

Water 

Hydrogen 

Oxygen 

6.677 xlO' 2 
3.339 xlO' 2 

6.677 XlO 22 
3.339 XlO 22 

Beryllium 


1.228 xlO" 1 

1.228 XlO 23 

Aluminum 


6.02 xlO* 2 

6.02 XlO 22 


Cross sections of the core calculational model are shown in figures 3 and 4. The 
atom densities used in the calculation are shown in table I. 


Neutron and Primary Gamma Source Generation 

The energy distribution of neutron source particles is based on a uranium 235 fis- 
sion spectrum. The rejection technique used to select neutron energies from a portion 
of the fission spectrum as approximated by Cranberg is described in reference 3. For 
primary gamma-energy selection, the reactor operating history shown in figure 5 was 
used to obtain the delayed gamma contribution using built-in delayed gamma tables. 
These gamma tables are described in reference 4. The spatial distribution of both neu- 
trons and primary gammas is determined based on the detailed axial and radial power 
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Local power density/average core power axial density 



TABLE H. - MEASURED RELATIVE POWER 
DENSITY AT MIDPLANE OF TWMR 
CRITICAL ASSEMBLY 


Fuel element 

Relative power density 

K-l 

0.744 

J-l 

.724 

J-2 

.712 

J-3 

.800 

H-l 

. 771 

H-2 

.671 

H-3 

.781 

H-4 

.883 

H-5 

.949 

G-2 

.709 

G-3 

.733 

G-4 

.851 

G-5 

.938 

G-6 

.983 

G-7 

1.000 


distributions (fig. 6 and table II) and fuel-element geometry. The angular distribution of 
the source particles was isotropic. Energy and region importance sampling were em- 
ployed in the calculation. This technique modifies the a priori energy and spatial distri- 
bution and, in general, can result in a reduced variance of the Monte Carlo answers for a 
given number of source particles. 


Secondary Gamma Source 

Secondary gamma source particles were generated at the sites of neutron inelastic 
scattering and absorption. These events were recorded on an interaction tape during the 
neutron Monte Carlo calculation. The information recorded included the type of inter- 
action, spatial location, energy of the neutron, and a number identifying the nuclide in- 
volved in the interaction. The interaction tape was then processed through the program, 
GASP, that generates the secondary gamma sources. The gamma production cross sec- 
tions that were used are discussed in reference 4. 

Monte Carlo Calculations 

Eighty thousand primary gammas were tracked to determine this contribution to the 
total gamma dose. Gammas produced from tracking 8000 neutrons were used to deter- 
mine the secondary gamma contribution to the total dose. A total of 144 tracking regions 


8 


were used in this calculation. The following table shows the time required on the 
IBM-7094 computer to perform the Monte Carlo Calculations: 


Type of source 

Computer time, 
min 

Primary gamma 

42.4 

Neutron 

51.3 

Secondary gamma 

41.3 


The cross section data used in ATHENA were specified in 81 equal lethargy intervals 
from 0.037 eV to 18.02 MeV for neutrons and from 10 keV to 10 MeV for photons. 


Results of Calculation 

The total dose calculated for the 16 chambers are shown in table III. The dose is 
based on the calculated energy deposition rates which were assumed to be constant over 

TABLE m. - COMPARISON OF CALCULATED AND MEASURED GAMMA DOSE 


Chamber 

location 

(a) 

, Distance 
from core 
bottom 
plate to 
center of 
chamber, 
cm 

Chamber 

Calculated gamma dose 
in graphite 

Measured gamma dose 
in graphite 

Percent 

deviation 

rads 

J/kg 

rads 

J/kg 

G-4 

84.7 

C-7 

680± 72 

6.8 ±0.72 

883 

8.83 

-23 

G-4 

62.1 

C-6 

1170± 94 

11.7 ± .94 

1253 

12.53 

+7 

G-4 

37.1 

C-l 

1153± 92 

11. 53+ .92 

1207 

12.07 

-4 

G-4 

9.5 

C-5 

97 5± 79 

9.75+ .79 

880 

8.80 

+11 

G-2 

62.1 

C-l 

987+104 

9. 87±1. 04 

839 

8.39 

+18 

G-2 

37.1 

C-ll 

744± 81 

7. 44+ .81 

881 

8.81 

-15 

K-l 

37. 1 

C-ll 

831± 97 

8. 31± .97 

850 

8.50 

-2 

g-10 

84.7 

C-13 

849± 89 

8. 49+ .89 

870 

8.70 

-2 

g-10 

37.1 

C-7 

1130± 93 

11.3 ± .93 

1279 

12.79 

-11 

g-10 

6.6 

C-10 

796± 86 

7. 96+ .86 

828 

8.28 

-4 

g-8 

37.1 

C-6 

1290+109 

12.9 ±1.09 

1201 

12.01 

+7 

g-6 

105.5 

C-13 

301+ 44 

3. 01± .44 

250 

2.50 

+20 

g-6 

62.1 

C-12 

960+ 78 

9.6 ± .78 

989 

9.89 

-3 

g-6 

37.1 

C-5 

1143+ 94 

11. 43± .94 

1055 

10.55 

+8 

g-6 

6.6 

C-3 

648± 74 

6. 48± .74 

723 

7.23 

+10 

g-4 

37.1 

C-1A 

87 5+ 84 

8. 75± .84 

804 

8.04 

+9 


a See fig. 10. 
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the 40-minute reactor operating time. A comparison between the calculated and the cor- 
responding experimentally measured values is given. The experimental measurements 
are discussed in the next section. 


GAMMA DOSE MEASUREMENTS 
Description of Graphite Ionization Chambers 

The graphite wall, carbon dioxide filled chambers were 0. 490 inch (1.24 cm) in 
diameter and 1.90 inch (4.83 cm) long and could be inserted within the control tubes or 
within the center post of the fuel element. The graphite chamber is shown in figure 7. 



1, lonizable volume - 0. 176 cm^ (carbon 

dioxide gas) 

2, Carbon electrode sheath (0.064 in. diam) 

3, Aluminum shield (0. 015 in. thick) 

4, Condenser (6000 (pF)) 

5, Electrode clip 

6, Aluminum condenser core 

7, Retainer ring 

8, Bellows assembly 

9, End cap 


10, Aluminum shell 1.900 in. long -490 in. 

diam - 0.015 in. wall 

11, Aluminum electrode cap 

12, Carbon chamber (0. 134 in. wall) 

13, Electrode insulator 

14, Aluminum electrode 

15, Electrode contact bar 

16, Condenser insulator 

17, Glass seal 

18, Contact pin 

19, End sleeve 


Figure 7. - 2500 rad graphite ionization chamber. 


These chambers were hermetically sealed and were charged and read by pushing the con- 
tact pin thereby stretching the bellows at the end of the chamber until electrical contact 
is made with the inner electrode. The electrical leakage of the chambers was about 
2 percent of full scale per day. 

Since the diameter of the chambers was restricted to permit insertion into the poison 
tubes, the areal density of the graphite walls was limited to 0. 544 gram per square cen- 
timeter. This thickness exceeds the range of a 1.25 MeV electron and, therefore, cor- 
responds to charged particle equilibrium for at least a 1. 5 MeV photon. Actually, 
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charged-particle equilibrium is generally achieved at a wall thickness considerably less 
than the range of the highest energy secondary electron produced by the primary photon 
interaction. The effective wall thickness is greater since the chambers are surrounded 
by aluminum when positioned in the fuel elements and by aluminum and water when posi- 
tioned within the poison tubes. If only aluminum surrounded the ionizable volume, then 
the measured dose would be 4 percent less than for carbon. If only water surrounded the 
ionizable volume, then the measured dose would be 11 percent greater. The largest 
error would therefore occur when the graphite wall chambers are surrounded by water. 

In this case, the total error caused by the walls not being in charged particle equilibrium 
for photons greater than 1.5 MeV should be no more than 4 percent. 

The graphite ionization chambers have a small inherent sensitivity to fast neutrons 
caused by carbon recoils. Since the carbon recoil atoms generated in the walls would be 
largely absorbed in the walls, the neutron response of the chamber would be largely that 
of the recoils in the carbon dioxide gas. The experimental results were analytically cor- 
rected for the neutron recoils in the carbon dioxide gas as discussed in the section Cor- 
rection for the Neutron Sensitivity of Graphite Chambers. 


TABLE IV. - LANDSVERK ELECTROMETER 
COMPANY CALIBRATION CERTIFICATE a 
[Range, 2500 rads(25J/kg) (carbon).] 


Serial 

num- 

ber 

Response of graphite dosimeters 

Cobalt 60 
gamma-ray 
source 

Gamma energy levels 

120 keV 

80 keV 

46 keV 

Cl 

76.2 

81.1 

89.6 

46.6 

C2 

78.4 

82.5 

90.3 

44.3 

C3 

94.4 

93.5 

88.8 

51.2 

C4 

74.0 

73.3 

66.9 

35.9 

C5 

86.0 

84.8 

89.6 

43.7 

C6 

84.4 

83.6 

77.7 

42.9 

C7 

87.0 

86.0 

81.6 

46.0 

C8 

84.8 

84.2 

76.8 

42.0 

C9 

80.0 

79.4 

82.7 

40.0 

CIO 

90.2 

89.3 

81.9 

43.5 

Cll 

85.6 

84.8 

80.3 

41.7 

C12 

85.2 

84.4 

77.3 

41.7 

C13 

107.2 

106.2 

96.0 

49.2 

C14 

85.8 

85.0 

81.4 

43.2 

C15 

84.4 

83.4 

75.5 

42.4 

C16 

94.0 

93.1 

85.4 

44.3 


a Dated Nov. 1965. 
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Figure 8. - Typical 2300 rad graphite ionization chamber response curve. 


The energy response for each graphite chamber obtained with a cobalt 60 source and 
with X-rays are given in table IV. Figure 8 shows a typical energy response curve for 
the graphite chamber. 


Absolute Calibration of Graphite Ionization Chambers 

The graphite chambers were calibrated at General Atomic by comparison with a 
standard thimble chamber in a bremsstrahlung beam resulting from 7 MeV Linac gener- 
ated electrons impinging on a thick fansteel (89 percent tungsten, 7 percent nickel* and 
'4 percent copper) target. The bremsstrahlung beam approximated a fission spectrum. 
Since all the photons were below the (y, n) threshold of the target elements, the beam 
contained no neutrons. 

The bremsstrahlung beam was directed to the center of a 4-inch- (10. 16 cm) thick, 
12-inch- (30.48 cm) diameter graphite disk shown in figure 9 containing eighteen 1/2- 

O 

inch- (1.27 cm) diameter, 2^ -inch- (6.985 cm) deep holes drilled on a 9 -inch (22.86 cm) 
diameter. All the graphite chambers were inserted into these holes and the graphite disk 
rotated at 2 revolutions per minute during the calibration. The thickness of the graphite 
between the ionization chamber and the front surface of the graphite was 1^ inches 
(3.175 cm) and thus slightly greater than the range of a 7 MeV electron; therefore, par- 
ticle equilibrium is established for photons at or below 7 MeV. 

All the graphite ionization chambers were intercalibrated on a relative basis during 
a single measurement. However, since the range of the graphite ionization chambers was 
several thousand roentgens and the standard calibration chamber was limited to 25 roent- 
gens (2.09 J/kg-air), a photodiode -plastic fluor detector (ref. 5) was used to monitor the 

3 

measurements during the intercalibration. The instantaneous dose rate was 1x10 roent- 
gens per second (8. 38 J kg"* (air) sec"*)) using a pulse length of 4. 5 microseconds and a 
pulse rate of 180 per second from the Linac. A 30-minute run resulted in about a half- 
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Figure 9. - Ionization chamber calibration apparatus. 

scale reading on the graphite chambers which corresponded to about 2600 roentgens (21.7 
J/kg-air). 

The standard calibration chamber was then inserted in an especially enlarged hole in 
the graphite disk and irradiated together with the photodiode plastic fluor detector until 
an approximate midscale reading on the standard calibration chamber was obtained. The 
standard calibration chamber behaved as an air -equivalent dosimeter from 20 keV to 
7 MeV. The exposure level of the standard chamber was compared with the reading on 
the photodiode -plastic fluor detector several times during this irradiation to obtain an 
average value. The maximum deviation between the several measurements of the aver- 
age was found to be less than 1 percent. 

The ratio of the photodiode plastic fluor detector reading obtained for the calibration 
of the graphite chamber to the reading obtained for the intercomparison with the standard 
chamber was multiplied by the standard chamber reading in roentgens (or J/kg-air) and 
then converted to rads (or J/kg) using the conversion factor 1 roentgen = 0.87 rad 
(0.00838 J/kg-air = 0.0087 J/kg-graphite). 

The resulting rad (J/kg) value was then divided by the decimal fraction of full-scale 
read on each of the graphite chambers during the intercomparison to obtain the full-scale 
reading in rads (J/kg) dose to graphite for each individual chamber. The results of the 
absolute calibration of the graphite chambers are tabulated in table V. 
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TABLE V . - DOSIMETER CALIBRATION USING 7 . 0 MeV BREMSSTRAHLUNG RADIATION 
(a) Graphite ionization chamber; gamma response coefficient, 0.9009 


Ionization 

chamber 

number 

Chamber reading, percent of full 
scale 

Value ob 
standard 

tained fror 
calibratio 
imber 

n Average corrected 
n chamber reading, 
full scale value 

Uncorrected 

reading 

Drift 

correction 

Corrected 

reading 

cha 

rads 

J/kg 

rads 

J/kg 

C-l 

56.4 

0 

56.4 

2300 

23.00*1 



C-l 

51.0 



51.0 

2030 

20. 30 J 

4039 

40.39 

C-2 

63.5 



63.5 

2300 

23.00 

3622 

36.22 

C-3 

63.5 



63.5 

2300 

23.00*1 



C-3 

57.6 



57.6 

2030 

20. 30 J 

3573 

35.73 

C-1A 

63.2 



63.2 

2300 

23.00 

3639 

36.39 

C-5 

57.9 



57.9 

2300 

23.00*1 



C-5 

51.9 



51.9 

2030 

20. 30J 

3941 

39.41 

C-6 

56.0 



56.0 

2300 

23.00 

4107 

41.07 

C-7 

58.0 



58.0 

2300 

23 . 00*1 



C-7 

52.2 



52.2 

2030 

20.30J 

3928 

39.28 

C-8 

63.2 



63.2 

2300 

23.00 

3639 

36.39 

C-9 

56.8 



56.8 

2300 

23.00*1 



C-9 

50.4 



50.4 

2030 

20. 30J 

4038 

40.38 

C-10 

61.4 



61.4 

2300 

23.00 

3746 

37.46 

C-ll 

56.2 



56.2 

2300 

23. 00*1 



C-ll 

49.7 



49.7 

2030 

20.30 j 1 

4089 

40.89 

C-12 

60.9 



60.9 

2300 

23.00 

3777 

37.77 

C-13 

71.5 



71.5 

2300 

23.00*1 



C-13 

64.0 



64.0 

2030 

20. 30 j 

3194 

31.94 

C-14 

58.9 



58.9 

2300 

23.0o1 



C-14 

52.5 



52.5 

2030 

20. 3 0J 

3886 

38.86 

C-15 

61.8 

3. 

7 

58.1 

2300 

23.00 

3959 

39.59 

C-16 

65.2 

.6 

64.6 

2300 

23.00 

3560 

35.60 


(b) Polyethylene ionization chamber; zero drift correction 


Ionization 

Chamber reading, percent 

Value obtained from 

Average corrected 

Gamma 

chamber 

number 

of full scale 

standard calibration 

chamber reading, 

response 

Uncorrected 

Corrected 



chamber 


full scale value 

coeffi- 


reading 

reading 

rads 

J/kg 

rads 

J/kg 

cient, 

€ 

P-1 

26.9 

26.9 

2630 

26 

30 

12 860 

128.60 

1.19 

P-1 

25.0 

25.0 

2320 

23.20 

12 860 

128.60 

1.25 

P-2 

22.2 

22.2 





14 920 

149.20 

1.29 

P-3 

22.3 

22.3 





14 090 

140.90 

1.22 

P-4 

23.6 

23.6 





14 190 

141.90 

1.30 

P-5 

22.0 

22.0 





14 220 

142.20 

1.21 

P-6 

25.5 

25.5 





11 400 

114.00 

a l. 13 

P-7 

26.8 

26.8 





12 150 

121.50 

1.26 

P-8 

28.0 

28.0 





12 380 

123.60 

1.35 

P-9 

22.0 

22.0 





14 620 

146.20 

a l. 64 

Average 









1.26 


^ot used in average. 



The possible individual errors in this calibration are 

(1) ±2 percent in the average energy to produce an ion pair in air for the roentgen to 

rad (J/kg-air to J/kg- graphite) conversion 

(2) ±2 percent in the reading of the graphite chambers 

(3) ±2 percent in the reading of the standard calibration chamber 

(4) ±1 percent in the value used to convert the standard calibration chamber reading 

to the intercalibration run for the graphite chambers. 

(5) ±2 percent in the photodiode -plastic fluor detector readings 

(8) ±2 percent in the absolute roentgen value of the standard calibration chamber 
It is estimated that this absolute calibration procedure gives the true absorbed 
gamma dose within ±5 percent. 


Positioning of Graphite Chambers in TWMR Critical Assembly 

As can be seen in figure 1, the TWMR critical assembly has a regular hexagonal 
cross section. Therefore, measurements in only one twelfth sector were required to ob- 
tain to an overall description of the gamma dose in the core. The positioning of the 
chambers is schematically shown in figure 10. The chambers were positioned above and 
below the zirconium stud using 3/8-inch- (0.952-cm) diameter, 0. 060-inch- (0.1524-cm) 
wall-thickness aluminum tube spacers. In the control tubes, 5/16-inch- (0.7938-em) 
diameter, 0.060-inch- (0.1524-cm) wall-thickness tube spacers were used and the space 
between the chambers filled with the cadmium nitrate solution. Also shown in figure 10 
are the locations of the polyethylene wall, ethylene filled ionization chambers used in 
connection with the total neutron plus gamma dose measurements described later. 

The measurements were made in four runs, each run employing 16 graphite wall and 
9 polyethylene -wall ionization chambers. For each run the same graphite-wall chamber 
was placed in the G-7 (center) fuel element and the same polyethylene -wall chamber in 
the g-16 control tube to monitor the relative power level between runs. The relative 
variation in the integrated dose and therefore the power level from run to run as indica- 
ted by the monitor ionization chambers was within 3 percent. 


Irradiation of Chambers in TWMR 

Each of the four irradiating runs lasted 40 minutes at full power. The reactor was 
brought to full power on a 30-second period. 

A measurement of the gamma time history was made for each of the four reactor 
runs using a gamma scintillation detector which discriminated against fast neutrons. 
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This detector was located about 2 feet (61 cm) above the beryllium reflector at the edge 
of the core and the relative gamma intensity measured as a function of time for each of 
the reactor runs . The relative gamma intensity as a function of time for the four runs 
are shown in figure 11. 





(d) Run 4. 

Figure 11. - Gamma intensity time history. 



Correction for Neutron Sensitivity of Graphite Chambers 

As has been noted before, the graphite chambers have an inherent sensitivity to fast 
neutrons because of the ionization caused by recoiling carbon and oxygen atoms in the 
carbon dioxide fill gas and, to a smaller extent, because of recoiling carbon atoms emit- 
ted from the graphite wall. The second effect is considered to be negligible since the re- 
coil atoms would be largely absorbed in the walls . 

The reading (in rads (J/kg)) of the graphite chamber R c resulting from irradiation 
in the TWMR critical assembly therefore consists of the dose KG in the graphite wall 
deposited by the photons plus the dose yN in the carbon dioxide gas deposited by the fast 
neutrons . 


R c = KG + i/N (1) 

(Symbols are defined in the appendix) where G is the gamma dose (rads or J/kg) depos- 
ited in water at the position of the graphite chamber, K is the ratio of gamma dose (rads 
or J/kg) in graphite to the gamma dose (rads or J/kg) in water. (The value of K is 
0.9009 (ref. 2)), N is the fast neutron dose (rads or J/kg) deposited to water at the posi- 
tion of the graphite chamber, and v is the ratio of the fast neutron dose in carbon dioxide 
to the fast neutron dose in water. The value of v was calculated as outlined later. 

When equation (1) is solved for the corrected gamma dose KG, and the ratio 
G/(G + N) is replaced by the equivalent spectrum index S, equation (2) results. 



This equation is used to correct the measured graphite chamber reading R^ to the 
true rad (J/kg) dose deposited in graphite KG. 

The ratio S was measured in 23 locations in the TWMR core and an average value 
for S of 0. 29±0. 03 obtained. The construction and calibration of the polyethylene cham- 
bers and the techniques used to obtain the ratio S for the TWMR core are discussed in 
the section TOTAL RADIATION DOSE MEASUREMENTS. 

The value of v was calculated from the following considerations. Fast neutrons 
colliding with low-mass atoms impart appreciable energy to the recoiling atom since the 
average fraction f of energy imparted per elastic collision (assuming isotropic scatter- 
ing in the center -of -mass system) is 2M/(M + l) 2 where M is the atomic mass of the 
recoil atom. The energy imparted to 1 gram of atoms D for a neutron flux per energy 
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initial AE centered about the energy E is 


D = ^or s fE(p(E)AE (3) 

where Jf is the atomic number density in atoms/gram, a g is the elastic scattering 
cross section, and <p(E)AE is the flux in the interval AE centered about the energy E. 
Thus, the total dose to 1 gram of material per unit time from fast neutrons is 



where E, is the lower limit of the fast neutron range. Since the chamber reading is 
proportional to the ionization caused by the recoiling nuclei, the lower limit of the fast 
neutron range is assumed to occur at the energy at which the recoiling atom moves so 
slowly that ionization of the material becomes improbable. The lower limit is believed to 
be 10 keV for the proton recoils (ref. 6) resulting from the bombardment of hydrogen by 
20-keV neutrons. The equivalent lower limits of the fast neutron range of 860 keV for 
carbon and 1. 5 MeV for oxygen are obtained by scaling up the hydrogen values consider- 
ing the lower average energy transfer and higher energy Required for equal speeds for 
these atoms. 

Thus, the expression for v, the ratio of the neutron dose per gram of carbon dioxide 
to the neutron dose per gram of water, is 


29 

z 
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(5) 


where C, O, and H are carbon, oxygen, and hydrogen, respectively. The integral of 
equation (4) has been replaced by a summation over the GAM II fine -group structure. The 
GAM II code (ref. 7) was used to supply fine-group elastic scattering cross section cr . 

Si 

and fine-group flux values <p. starting at 14.9 MeV for the calculation. The AU. are 
the lethargy widths for the fine groups (AU = 0.1 from 14.9 MeV to 86. 5 keV; AU = 0.25 
from 86. 5 keV to E^), and E. is the energy at the midpoint of each interval. 

The shape of the neutron spectrum <p (U) computed by the GAM II code is dependent on 
the input buckling. The fast group buckling obtained from a one -dimensional transport 
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theory TDSN (ref. 8) calculation of the TWMR critical assembly was used in GAM II. The 
buckling dependent, zero dimensional GAM n computed neutron spectrum is compared 
with the ATHENA computed three-dimensional neutron spectrum in the center water re- 
gion in figure 12. Good agreement is noted both in the spectrum and in the magnitude of 
the flux. 



Figure 12. - Comparison of GAM II and Athena calculated neutron spectrum in central water region of TWMR 
critical aseembly. 


The value computed for v using equation (5) was 0. 0634±0. 0187. The error assigned 
to v reflects the lack of reliable data on the lower energy limit of the fast neutron range 
for hydrogen, carbon, and oxygen and to a much lesser extent (because the flux appears 
in both the numerator and denominator of eq. (5)) on the uncertainty in buckling value 
used in the GAM n spectrum calculation. The magnitude of the error was obtained by re- 
calculating v using a GAM II spectrum in which the buckling was reduced by 30 percent 
to the overall epithermal buckling value used in neutronic calculations for this core. The 
lower limit of the fast neutron range for all the elements was taken to be 19.3 keV. The 
buckling change has only a few percent effect on v, while the reduced lower limit of the 
fast neutron range tends to greatly increase v. The recalculated value for v was 0. 1008. 
This value is considered to be the upper limit of the ratio of neutron dose in carbon diox- 
ide to the neutron dose in water. Since the uncertainty is one-sided about v, one-half of 
this difference in v was assigned as the error in v. It must be noted that although the 
error assigned to v is nearly 30 percent, the resulting error in KG, the gamma dose in 
graphite as obtained by equation (2), is only 4 percent. 
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Results of Gamma Dose Measurements 


The absolute gamma dose in the graphite KG was obtained from the graphite cham- 
ber readings Rj-, using equation (2) . Substitution of the parameters u, K, and S leads 
to equation (6). 


KG = 


1 

= R C 

1 

1 +-(=-') 

x ( 0.06341 / 1 A 

K \S /J 


0.9009 \0. 29 /_ 


= 0. 853 Rq 


( 6 ) 


The data and results are listed in tables HI and VI. 


Error Analysis of Gamma Dose Measurements 

The relation between the parameters in equation (2) was used to propagate the un- 
certainties in v, S, and R c to determine the resulting error in the measured gamma 
dose to graphite. 



The uncertainty in R^ consists of the ±5 percent uncertainty in the calibration of the 
graphite chambers discussed earlier, a ±2 percent uncertainty in reading the graphite 
chamber, a ±3 percent uncertainty in the time at which the reactor was at power, and a 
±7 percent uncertainty in the measurement power level of the reactor which is discussed 
in the section Measurement of Absolute Power Level. Thus, the total uncertainty in R^-, 
is ±9 percent. If the error in measurement of the power level is neglected, the uncer- 
tainty is ±6 percent. As was discussed earlier, an uncertainty of ±30 percent was as- 
signed to v, and an uncertainty of ±10 percent resulted from the measurements of S in 
the TWMR. Thus, the total error in the value of KG as given by equation (7) is ±10 per- 
cent. 
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TABLE VI. - RESULTS OF GRAPHITE IONIZATION CHAMBER MEASUREMENTS IN 

TWMR CRITICAL ASSEMBLY 


Control tube 

Distance from 

Chamber 

Chamber reading with run 

Absorbed gamma dose 

or fuel 

core bottom 

number 

normalization included, 

in graphite, 

element 

plate to center 



R c 

KG a 


of chamber, 








cm 



rads 

JAg 

rads 

JAg 

G-2 


.4 

C- 

16 

698 

6.98 

594 

5.94 


37.1 

C- 

11 

1035 

10.35 

881 

8.81 


62.1 

c- 

1 

986 

9.86 

839 

8.39 


105.1 

c- 

12 

362 

3.62 

309 

3.09 

G-4 

9.5 

c- 

5 

1034 

10.34 

880 

8.80 


37.1 

c- 

1 

1418 

14.18 

1207 

12.07 


62.1 

c- 

6 

1472 

14.72 

1253 

12.53 


84.7 

c- 

7 

1038 

10.38 

883 

8.83 


105.1 

b c- 

15 

194 

1.94 

165 

1.65 

G-7 

9 

.4 

c- 

5 

1213 

12.13 

1032 

10.32 


37 

.1 

c- 

12 

1717 

17.17 

1461 

14.61 


62 

.1 


r Run 1 

1675 

16.75 

1426 

14.26 





Hun 2 

1679 

16.79 

1429 

14.29 




C-14-j 

Run 3 

1678 

16.78 

1428 

14.28 





Run 4 

1678 

16.78 

1428 

14.28 


84.7 

b c- 

15 

930 

9.30 

792 

7.92 


105 

.1 

b c- 

15 

202 

2.02 

172 

1.72 

J-3 

9 

.4 

c- 

13 

1034 

10.34 

880 

8.80 


62 

.1 

c- 

12 

1326 

13.26 

1128 

11.28 


105 

.1 

c- 

16 

402 

4.02 

342 

3.42 

K-l 

9.4 

c- 

16 

614 

6.14 

522 

5.22 


37 

.1 

c- 

11 

999 

9.99 

850 

8.50 


62 

.1 

c- 

11 

965 

9.65 

821 

8.21 


84.7 

c- 

1 

709 

7.09 

603 

6.03 


105 

.1 

c- 

13 

338 

3.38 

288 

2.88 

g“4 

6 

.6 

c- 

3 

792 

7.92 

674 

6.74 


37 

.1 

c- 

1A 

945 

9.45 

804 

8.04 


62 

.1 

c- 

1 

954 

9.54 

812 

8.12 


110 

.0 

c- 

3 

236 

2.36 

201 

2.01 

g-6 

6 

.6 

c- 

3 

850 

8.50 

723 

7.23 


37 

.1 

c- 

5 

1239 

12.39 

1055 

10.55 


62.1 

c- 

12 

1162 

11.62 

989 

9.89 


84 

.7 

b C- 

15 

543 

5.43 

462 

4.62 


110 

.0 

c- 

13 

180 

00 

o 

154 

1.54 


a KG = 0.853 R c . 

^Suspected defective ionisation chamber. 



TABLE VI. - Concluded. RESULTS OF GRAPHITE IONIZATION CHAMBER 
MEASUREMENTS IN TWMR CRITICAL ASSEMBLY 


Control tube Distance from Chamber Chamber reading with run Absorbed gamma dose 


or fuel 
element 


core bottom 
plate to center 
of chamber, 


number normalization included, 
Ro 


in graphite, 



*KG = 0.853 R c . 

^Suspected defective ionization chamber. 


COMPARISON OF CALCULATION AND MEASUREMENT OF ABSOLUTE GAMMA DOSE 


The measured and calculated absolute gamma dose deposited in graphite are com- 
pared in figure 13 and in table HI. For convenience in comparing the calculated and 
measured values, the solid lines shown in figure 13 enclose the band of probable error 
of the measured values. The calculated values are shown as bars indicating the one 
standard deviation confidence level. It is noted that the measured and calculated values 
agree within their mutual uncertainties in all but one case. Table HI shows that the de- 
viation between experiment and calculation is within 10 percent in 10 out of 16 locations 
compared. 


■g 



Distance from core bottom plate, cm 
(b) Axial profile fuel element, G-4. 

Figure 13. - Comparison of measured and calculated absolute gamma dose. 
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Athena Monte Carlo calculated 
Probable error limit measured 

i i i i l_ 

(c) Axial profile control tube, g-6. 



J I i i i 

30 40 50 60 70 

Distance from core bottom plate, cm 

(d) Axial profile control tube, g -10. 

Figure 13. - Continued. 





Radial distance from core center, cm 
(f) Radial profile in fuel assemblies 37. 1 centimeters above bottom plate. 
Figure 13. - Concluded. 


NEUTRON DOSE CALCULATIONS 

The neutron dose in water is the result of elastic collisions of neutrons with the 
nuclei in the water. The recoiling nuclei deposit their kinetic energy as heat. The other 
processes by which neutron heating can occur are inelastic scattering, neutron capture, 
and charged particle reactions. Inelastic scattering is nonexistent in hydrogen but has a 
threshold above 6 MeV in oxygen. The inelastic scattering gamma rays were included in 
the gamma calculation but the small amount of kinetic energy deposited by the recoiling 
neutrons was ignored in the neutron dose calculations. The neutron capture gammas 
were included in the gamma calculations and the charged particle reactions were assumed 
to be negligible. 

The method used to calculate neutron dose was outlined in connection with the calcu- 
lation of v in the section Correction for Neutron Sensitivity of Graphite Chambers. The 
total dose D tot (r ) or energy deposited in 1 gram of material per unit time from fast neu- 
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trons at the position r is 



Jfa fE<p(E, r )dE 
s 


( 8 ) 


where is the atomic number density in atoms/gram, o is the elastic scattering 
cross section of the atom, f is equal to 2M/(M + 1) where M is the atomic mass, and 
<p(E, r ) is the neutron flux per unit energy centered about the energy E at the position of 
measurement r with the reactor operating at power. 

In the calculation, the integral of equation (8) is replaced by a summation over a 
30 group GAM n structure. The upper limit on the summation was taken at 14. 9 MeV. 
The lower limit E^ was taken at 854 eV. The contribution to from neutrons out- 

side this energy range is negligible. 


HoO _ 

D t0 2 t (r)-A(r) 


j =30 




i=i 


( 9 ) 


The GAM II code is used to supply the group elastic scattering cross sections a .. 

S J 

The flux values cp ^ are obtained from a radial three -region P Q ®4 transport unit cell 
calculation using the TDSN program (ref. 8). The central radial region in this calculation 
consisted of the homogenized uranium 235, uranium 238, aluminum, and tungsten in a fuel 
element. The next radial region contained the homogenized aluminum of the pressure tube 
and the two control tubes associated with the fuel element, and the cadmium within the 
control tubes. The outer radial region in the calculation contained the water in the mod- 
erator and in the control tubes associated with the fuel element. Thirty group spectrum 
averaged cross sections supplied by the GAM n code were used in each region. The 
TDSN calculated fluxes for each region are normalized to a source strength of 1 neutron 
per square centimeter of the cell volume. The multiplier A(r) is used to renormalize 
the flux to that in the TWMR at the position r at power. 

As previously stated, the shape of the GAM II computed neutron spectrum is depen- 
dent on the value of the input buckling. The input buckling used in the calculation of v 
was used in this calculation. 

The value calculated for the quantity within the summation sign in equation (9) was 

4.399 Me Y ^51^ = 7.048X10 -10 J_cm^ 

g(H 2 0) fission neutron kg(H 2 0) fission neutron 
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for the water region in the TDSN radial cell. Similarly, a value of 


8.053X10' 10 


J cm^ 

kgfHgO) fission neutron 


was calculated for the fuel region. 

The normalizing factor A(r ) was computed using the power distributions measured 
in the TWMR critical assembly. These power distributions were in good agreement with 
neutronic calculations. The power distributions in the TWMR critical assembly were 
separable; that is, the shape of the power distribution in the axial direction was indepen- 
dent of the radial position of the fuel element along which the axial measurements were 
made. The axial power shape used in the calculation is shown in figure 6. The radial 
power density at the midplane of the core normalized to unity at center fuel element is 
shown in table n. 

The value A(r) is the product of a radial factor F^(r) corresponding to the relative 
power density at the radius r at which the neutron dose is to be calculated, and an axial 
factor F^(z) corresponding to the power density at a distance z above the core bottom 
plate at which the calculation is made, and a constant factor corresponding to the power 
level of the critical assembly during the measurement. 

The radial factor Fp^(r) when the ionization chamber is located within the fuel ele- 
ment is taken directly from table II. If the ionization chamber is located in a control 
tube, then the radial factor is the average value of the three adjacent fuel elements. The 
axial factor F^(z) is found using figure 6. 

The constant factor is calculated using the absolutely measured total core power level 
of 84.2 watts and the 40 minute irradiation time and the 657 100 cubic centimeter TWMR 
core volume as follows: 

84. 2 W x 3. 27X10 10 ^- s - si o n x 2. 477 - so urce neutrons x 4Q min x 60 se c 

W-sec fission min 

x 1 _ 2 491x10*® source neutrons 

657 100 cc cc 

Thus, A(r) is 2.49x10*® F^(r) x F^(z) and equation (9) becomes, in the water 
regions, 
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I 


D t ot (r) = 1.756X10 1 F R (r) x F A (z) J 


kg(H 2 0) 

= 1.756X10 3 F R (r) x F A (z) (rads) (10) 

Similarly, in the fuel regions 

D tot (r) = 2.006X10 3 F R (r) X F A (z) rads; 2.006X10 1 F R (r) x F A (z) J/kg(H 2 0) 


RESULTS AND ERROR ANALYSIS OF NEUTRON DOSE CALCULATIONS 

The calculated neutron dose is given in table VII. The errors in these calculations 
result from the uncertainty in the GAM II computed neutron spectrum caused by sensi- 
tivity of the spectrum to the buckling and the uncertainty in the power distribution used in 
the calculations. 

If the input buckling to the GAM II spectrum calculation is decreased by 30 percent to 
the overall epithermal buckling value used for neutronic calculations for the core, the 
value computed for the sum in equation (9) is increased by 7 percent. The uncertainty in 
the power distributions is about ±5 percent. The total error associated with the calcula- 
tion is therefore ±6 percent assuming that the error connected with the buckling input to 
GAM II is one-sided and, therefore, contributes only ±3. 5 percent to the total error. 


TOTAL RADIATION DOSE MEASUREMENTS 

The total radiation dose measurements were made using the gamma sensitive graphite 
chambers previously described and the neutron sensitive polyethylene chambers to be de- 
scribed. Both the graphite and the polyethylene chambers were absolutely calibrated in 
the bremsstrahlung gamma beam of the Electron Linear Accelerator (Linac) at General 
Atomic and gamma response coefficients (chamber dose reading per unit rad (J/kg) 
gamma dose in water) obtained. This calibration (for the graphite chambers) has already 
been described. 

The graphite and polyethylene chambers were then irradiated in the mixed neutron 
and gamma spectrum in the TRIGA reactor at General Atomic. The total radiation de- 
position to water was absolutely measured at the test location with the TRIGA reactor 
using a water-filled calorimeter. 
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TABLE VH. - RESULTS OF POLYETHYLENE IONIZATION CHAMBER CALCULATIONS 


AND MEASUREMENTS IN TWMR CRITICAL ASSEMBLY 


Control tube 
or fuel 
element 

Distance from 
core bottom 
plate to center 
of chamber, 
cm 

Calculated neutron dose 
in water 

Measured neutron dose 
in water 

Percent 

difference 

rads 

JAg 

rads 

JAg 


G-2 

9.4 

1760 

17.6 

— 



— 


37.1 

2260 

22.6 

— 

— 

— 


62.1 

2160 

21.6 

— 

— 

— 


67.0 

2060 

20.6 

2180 

21.8 

-5 


105.1 

739 

7.39 

— 

— 

— 

G-4 

9.5 

2120 

21.2 

— 

— 

— 


a 32. 4 

2650 

26.5 

2230 

22.3 

+19 


37.1 

2710 

27.1 

— 

— 

— 


62.1 

2590 

25.9 

— 

— 

— 


67.0 

2470 

24.7 

3660 

36.6 

-32 


84.7 

1840 

18.4 

— 

— 

— 


105.1 

888 

8.88 

— 

— 

— - 

g-4 

6.6 

1630 

16.3 

— 

— 

— 


11.4 

1530 

15.3 

1560 

15.6 

-2 


32.4 

1910 

19.1 

2390 

23.9 

-20 


37.1 

1970 

19.7 

— 

— 

— 


62.1 

1880 

18.8 

— 

— 

— 


67.0 

1790 

17.9 

2280 

22.8 

-22 


80.0 

1600 

16.0 

2000 

20.0 

-20 


105.3 

630 

63 

810 

8.10 

-22 


110.0 

444 

4.44 

— 

— 

— 

g-8 

6.6 

2060 

20.6 

— 


— 


11.4 

1930 

19.3 

2350 

23.50 

-18 


32.4 

2420 

24.2 

3200 

32.0 

-24 


37.1 

2490 

24.9 

— 

— 

— 


62.1 

2380 

23.8 

— 

— 

— 


80.0 

1880 

18.8 

2520 

25.2 

-25 


84.7 

1690 

16.9 

— 

— 

— 


105.3 

797 

7.97 

1130 

11.3 

-29 


110.0 

562 

5.62 

— 

— 


g-12 

6.6 

2290 

22.9 

— 

— 




a ll. 4 

2160 

21.6 

2060 

20.6 

+5 


a 32. 4 

2700 

27.0 

3410 

34.1 

-21 


37.1 

2760 

27.6 

— 

— 

— 


62.1 

2630 

26.3 

— 

— 

— 


a 67. 0 

2520 

25.2 

1620 

16.2 

+55 


a 80. 0 

2080 

20.8 

2970 

29.7 

-30 


84.7 

1870 

18.7 

— 

— 

— 


a 105. 3 

840 

84.0 

970 

9.7 

-9 


110.0 

630 

63.0 

— 

— 

— 


a Suspected defective ionization chamber. 



TABLE VH. - Concluded. RESULTS OF POLYETHYLENE IONIZATION CHAMBER 
CALCULATIONS AND MEASUREMENTS IN TWMR CRITICAL ASSEMBLY 


Control tube 
or fuel 
element 

Distance from 
core bottom 
plate to center 
of chamber, 
cm 

Calculated neutron dose 
in water 

Measured neutron dose 
in water 

Percent 

difference 

rads 

J/kg 

rads 

J/kg 

h-8 

6.6 

2090 

20.9 

— 

— 

— 


a ll. 4 

1970 

19.7 

1110 

11.1 

+78 


32.4 

2460 

24.6 

2960 

29.6 

-17 


37.1 

2530 

25.3 

— 

— 

— 


62.1 

2410 

24.1 

— 

— 

--- 


a 67. 0 

2310 

23.1 

1500 

15.0 

+53 


80.0 

1900 

19.0 

2590 

25.9 

-26 


84.7 

1710 

17.1 

— 

— 

— 


105.3 

810 

8.1 

860 

8.6 

+6 


110.0 

570 

5.7 

— 

— 

— 

j-4 

6.6 

1710 

17.1 

— 

— 

— 


11.4 

1620 

16.2 

1550 

15.5 

+5 


32.4 

2020 

20.2 

2490 

24.9 

-19 


37.1 

2070 

20.7 

— 

— 

— 


62.1 

1980 

19.8 

— 

— 

— 


67.0 

1880 

18.8 

2630 

26.3 

-28 


80.0 

1570 

15.7 

2000 

20.0 

-22 


84.7 

1410 

14.1 

— 

— 

— 


105.3 

663 

6.63 

1050 

10.5 

-37 


110.0 

468 

4.68 

— 

— 

— 

k-1 

6.6 

1670 

16.7 

— 

— 

— 


11.4 

1580 

15.8 

1570 

15.7 

0 


32.4 

1970 

19.7 

2270 

22.7 

-13 


37.1 

2020 

20.2 

— 

— 

— 


62.1 

1930 

19.3 

— 

— 

— 


67,0 

1840 

18.4 

2130 

21.3 

-14 


84.7 

1380 

13.8 

— 

— 

--- 


105.3 

647 

6.47 

724 

7.24 

-11 


110.0 

457 

4.57 

— 

— 

— 


a Suspected defective ionization chamber. 



The Linac and the TRIGA measurements were then combined with equation (1) to de- 
termine a gamma response coefficient e and a neutron response coefficient 6 for the 
polyethylene chamber analogous to the gamma response coefficient K and the neutron 
coefficient v for the graphite chambers given in equation (1). The gamma or neutron re- 
sponse coefficients are defined as chamber dose reading per unit gamma or neutron rad 
(j/kg) dose in water. 

Having thus obtained neutron and gamma response coefficients for the graphite and 
polyethylene chambers, it was possible to solve independently for the neutron dose to 
water at the calibration location in the TRIGA reactor, and at each measurement location 
in the TWMR reactor. The average ratio of gamma dose to total dose S for each reactor 
can thus be determined. 


Description of Polyethylene Ionization Chambers 


Nine polyethylene wall, ethylene -filled chambers were used. These chambers were 
sensitive to both gamma rays and fast neutrons. The chambers were 1. 85 inches 
(4. 67 cm) long and had an ionizable volume of 0. 040 cubic centimeter. The outer diam- 
eter of the chambers was 0. 490 inch (1.245 cm) for the same reasons as the graphite wall 
dosimeters. Figure 14 shows a schematic of the polyethylene chamber. The resulting 



1, Aluminum shell (L 850 in. long - 0. 490 in. 

diam -0.015 in. wall) 

2, Aluminum electrode (0.025 in. diam) 

3, Bakelite shield 

4, Condenser (3000 pF) 

5, Electrode clip 

6, Aluminum condenser core 

7, Retainer ring 

8, Bellows assembly 

9, End cap 


10, Ionizable volume 0.040 cm^ 

(ethyelene gas) 

11, Polyethelyene chamber 
(0. 165 in. wall) 

12, Electrode insulator 

13, Electrode contact bar 

14, Condenser insulator 

15, Glass seal 

16, Contact pin 

17, End sleeve 


Figure 14. - 5000 rad polyethylene ionization chamber details. 
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polyethylene wall thickness of 0. 385 gram per square centimeter exceeded the range of 
1.0 MeV and thus provided an equilibrium thickness for 1.25-MeV gamma rays. 

Since the polyethylene chambers were placed mainly in poison tubes, they were, in 
most measurements, surrounded by water which increased the effective wall thickness. 
The absorbed dose in water and polyethylene varies by only 3 percent in the Compton 
region and therefore the maximum error caused by the walls not being infinite was about 
3 percent. 


Calibration of Polyethylene Chambers 


First, the gamma response of the polyethylene chambers was absolutely determined 
using the bremsstrahlung beam of the Linac. This calibration was identical with that for 
the graphite chambers discussed above. The results of this calibration is given in 
table V and figure 15. 


120 r 


100 [ 


80 1 


.01 



I I I I 
1.0 

Photon energy, MeV 


j i i_i 


10 


Figure 15. - Typical polyethylene ionization chamber response curve. 


Second, the response of the polyethylene chambers to the mixed gamma plus neutron 
radiation spectrum within the General Atomic TRIGA reactor was absolutely measured by 
comparison against a specially designed water-filled calorimeter. Ideally, this calibra- 
tion should have been made in the mixed radiation spectrum of the TWMR core since the 
calibration would then be directly applicable to the data measured there. However, the 
allowable flux levels in the TWMR were insufficient to cause an appreciable rise in the 
water bath temperature of calorimeter, and so the calibration was made in the high- 
power density TRIGA core. It was therefore necessary to determine the ratio of the 
gamma to total gamma plus neutron dose S for each of the cores and apply suitable cor- 
rections to the data. 

Houghton, Jupiter, and Trimble (ref. 2) designed the calorimeter and made the dose 
measurements. The outer jacket of the calorimeter was held to a maximum diameter of 
1.250 inches (3. 17 cm) to allow it to fit into a tube in the TRIGA core. Since the size was 
restricted, they used an adiabatic jacket to reduce heat transfer between the water ab- 
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sorbing mass and its environment. This adiabatic jacket was placed halfway between the 
outer jacket and the water mass and was made of aluminum and epoxy with a heating coil 
potted into the epoxy; adjustment of the power dissipated in the coil almost completely 
eliminated heat transfer between the water mass and its environment. They further re- 
duced heat transfer by using fine (0.006 in. (0.0152 cm) diameter) copper connection 
wires to lead into the water mass and keeping all heat paths between the water mass and 
the adiabatic jacket long and of small cross sectional area. All interior surfaces were 
painted white and the area around the jacket was evacuated. Since the neutron and gamma 
energy would be absorbed by the vessel walls as well as by the water, they designed the 
vessel to have a small mass compared with the water mass. Sheet polyvinyl chloride 
(PVC) was thermoformed into a two-piece bottle with approximately 0.004-inch (0.010 cm) 
thick walls. A thin (~0. 002 in. (0. 0051 cm)) coating of epoxy was painted on the outside 
of the vessel to eliminate vapor pumping since the PVC is not impervious to water vapor. 
Details of the calorimeter are shown in figure 16. 



Adiabatic 

jacket (a) Calorimeter electronics block diagram. 


Figure 16. - Calorimeter details. 


Fine control of the temperature difference between the water mass and the adiabatic 
jacket and the ability to accurately sense small incremental temperature changes in the 
water mass dictated the use of thermistors as temperature sensing elements. The 
thermistors had a nominal resistance of 25° C of 50 000 ohms, a spherical diameter of 
0. 043 inch (0. 109 cm) and a temperature coefficient of resistance of 4. 5 percent per Cen- 
tigrade degree. A pair of these thermistors, matched to within 0. 2 percent of each other, 
was imbedded, one in the adiabatic jacket and one in the water mass: it forms two legs 
of a bridge circuit as shown in figure 16(c) and is used to control the temperature of the 
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Number 24 stranded leads- 


Copper feed-thrus, 
1/8 in. diam _ 


Poured epoxy top seal -% 
Lucite spacer- 

Soldered connection-^ 

Water mass copper lead 
wires, 0.006 in. diam_— ] 
0.010 in. polyvinyl- 
chloride (PVC) spacer-* 


0.005 in. Nichrome heater^ 
wire cast into jacket 

Water mass thermistor - 

0. 005 in. Nichrome 
water mass heater coil — "1 


0.020 in. PVC spacer - 


0.020 in. PVC spacer - 


To vacuum pump 


5/16 in. o. d. by 0.03 in. wall 
stainless steel tube 


Epoxy closure seal 

0.010 in. diam copper 
jacket lead wires 

0.210 in. o. d. by 0.037 in. 
wall PVC neck tube 

^ Jacket thermistor 


-Water level 


Poured epoxy bottom seal 

(b) Cross section of calorimeter. 
Figure 16(b). - Continued. 


1.375 in. 


0.577 in. o.d. PVC 
bottle walls, 0.004 in. thick with 
0. 002 epoxy outer coating 

|^0. 780 in. i.d. by 0. 20 in. wall aluminum 
jacket liner, 0.020 in. paper/epoxy coat- 
ing on outside of jacket 

‘M. 250 in. o.d. by 0.062 in. 
wall aluminum outer tube 
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adiabatic jacket with respect to the water mass to eliminate heat transfer. An additional 
thermistor in the water mass forms one leg of the other bridge circuit of figure 16. This 
was carefully calibrated to allow accurate monitoring of the rate of temperature rise in 
‘the water mass. Houghton and his coworkers (ref. 2) calibrated the thermistor as fol- 
lows: The entire water mass assembly was placed in a large volume of water contained 
by a vacuum -jacketed glass-walled flask. A copper -constantan thermocouple was located 
on the surface of the calorimeter water mass and, using an ice bath cold junction, con- 
nected to a precision potentiometer which allowed accurate determination of the tempera- 
ture of the water bath. A battery -driven resistance heater was used to raise slowly the 
temperature of the water bath and the calorimeter water mass assembly suspended within 
it. The output of the bridge circuit as a function of temperature change was monitored by 
a microvoltmeter and a strip chart recorder, and was found to be 14.85 millivolts per °C 
which agreed very well with calculations based on the manufacturer's rated temperature 
coefficient of resistance. This constant was used in the analysis of subsequent measure- 
ments. A 4.25-ohm heater coil was included in the water absorbing mass to control the 
temperature and provide a check on the operation of the system. 

Once the jacket controls were correctly adjusted, they (ref. 2) tested the calorimeter 
by applying power to the water mass coil and observing the rate of temperature rise of the 
.water mass on the strip-chart recorder. The slope of the temperature rise curve (in 
mV/min) gives the heating rate since the thermistor constant is known to be 14.85 milli- 
volts per °C. 

When the calorimeter was operated in the TRIGA reactor, a temperature rise rate of 
0.0829° C per minute was observed at a steady-state reactor power level of 10 kilowatts; 
this rate was obtained from the recorder trace which showed a voltage change rate (due to 
the change in resistance in the calibrated thermistor in the water mass) of 1.23 millivolts 
per minute. Since 0.0829° C per minute corresponds to 0.0829 gram-calorie per gram 
per minute for pure water at 13° C (the temperature of the reactor water during the cali- 
bration runs) the heat input rate to the calorimeter water mass is given by 

0 . 0829 g ~ cal 4. 19X10 7 ergs = 3. 47X10 6 ergs 
g-min g-cal g-min 


and the dose rate in water at 10 kilowatts is 


fi er S s 

3.47x10 _ a rads 9 * 

g ~ min = 3. 47X10* “TT (3. 47X10 Z J kg" 1 min" 1 ) 

10 2 ergs 
g-rad 


The reactor power level during the irradiation of the dosimeters was 1.8 kilowatts; 
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therefore the dose rate in water corresponding to the readings of the dosimeters is 6250 
rads per minute (62. 5 J kg”* min"*). 

As a check on the calorimeter operation, the reactor was operated at 20 kilowatts 
with a heating rate of 0. 167° C per minute being obtained in the calorimeter water mass. 
Doubling the power doubles the heating rate which indicates that the relative power rates 
of the reactor are well known and that the calorimeter is able to accurately follow the re- 
actor power level changes. 

The results of the absolute calibration of the polyethylene chambers against a water 
calorimeter are shown in table VUI. Also shown is the calibrated rad reading of several 
graphite chambers which were placed in the same reactor core position as the polyethyl- 
ene chambers. 


TABLE Vm. - CALIBRATION OF ION CHAMBERS IN TRIGA REACTOR 


(a) Polyethylene ionization chamber; TRIGA calorimeter dose 
6250 rads (62.5 J/kg) 


Ionization 

chamber 

number 

Chamber reading, percent of 
full scale 

Full-scale average 
corrected chamber 
reading 

Uncorrected 

reading 

Position 

correction 

Corrected 

reading 

rads 

J/kg 

P-1 

49.5 

0.982 

48.6 

12 860 

128.6 

P-2 

41. 1 

1.02 

41.9 

14 920 

149.2 

P-3 

43.4 

1.02 

44.3 

14 090 

140.9 

P-4 

45.0 

.982 

44. 1 

14 190 

141.9 

P-5 

43.1 

1.02 

43.9 

14 220 

142.2 

P-6 

55.4 

.982 

54.4 

11 400 

114.0 

P-6 

56. 3 

.982 

55.3 

11 310 

113.1 

P-7 

50.4 

1.02 

51.4 

12 150 

121.5 

P-8 

51. 5 

.982 

50.6 

12 360 

123.6 

P-9 

43.5 

.982 

42.7 

14 620 

146.2 

P-9 

41.3 

1.02 

42.2 

14 810 

148.1 

P-9 

42.4 

1.02 

43.3 

14 420 

144.2 


(b) Graphite ionization chamber 


Ionization 

Chamber reading, percent of 

Graphite chamber 

chamber 

full scale 


reading in TRIGA 

number 





if i 

Uncorrected 

Position 

Corrected 

reactor, Kq 


reading 

correction 

reading 

rads 

JAg 

C-2 

60.8 

1.02 

62.0 

2 246 

22.46 

C-5 

57.4 

1.02 

58.6 

2 309 

23.09 

C-9 

53.3 

1.02 

54.3 

2 193 

21.93 

C-14 

56.5 

.982 

55.5 

2 157 

21.55 

Average 

2 226 

22.26 
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Positioning and Irradiation of Polyethylene Chambers in 
TWMR Critical Assembly 


The positioning and irradiation of the polyethylene chambers in the TWMR critical 
assembly have been described in connection with the graphite chambers in previous sec- 
tions Positioning of Graphite Chambers in TWMR Critical Assembly and Irradiation of 
Chambers in TWMR. 


Analysis of Total Dose Data 

The gamma sensitive graphite chambers and the polyethylene neutron sensitive cham- 
bers responded to some extent to both neutrons and photons. Therefore, the total dose 
calibration of the chambers against the water -filled calorimeter in the mixed neutron and 
gamma radiation spectrum of the TRIGA reactor could not be directly applied to. the total 
dose measurements in the TWMR, which had a different mixed radiation field. The data 
measured in the Linac and the TRIGA reactor was therefore used to determine individual 
response coefficients for neutrons and for gammas for the graphite chambers and the 
polyethylene chambers. These response coefficients were then applied to the TRIGA and 
the TWMR critical assembly dose measurements to obtain the gamma dose and the neu- 
tron dose separately at each position measured. In addition average spectral indices, 
that is, the average ratio of the gamma dose in water to the total mixed gamma plus neu- 
tron radiation dose in water, were determined for the TRIGA reactor and for the TWMR 
critical assembly and were found to differ by 20 percent. 

Determination of chamber response coefficients . - The chamber response coeffi- 
cients e, 6, K, v are defined as follows: 

For the polyethylene chamber, 


Rp = eG + 6N 


(ID 


For the graphite chamber, 


R c = KG + i/N (12) 

where Rp is the reading of the polyethylene chamber for mixed radiation dose of G rads 
(J/kg) gamma dose plus N rads (J/kg) neutron dose measured in water, e is the dose 
reading of the polyethylene chamber per rad (J/kg) of gamma dose in water, G is the rad 
(J/kg) gamma dose at the position of measurement in water, 6 is the dose reading of the 
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polyethylene chamber per rad (J/kg) neutron dose in water, N is the rads (J/kg) neutron 
dose measured in water, R^, is the reading of the graphite chamber for the mixed radia- 
tion dose of G rads (J/kg) gamma dose plus N rads (J/kg) neutron dose measured in 
water, K is the response of the graphite chamber per rad (J/kg) gamma dose in water, 
and v is the response of the graphite chamber per rad (J/kg) neutron dose in water. 

In the Linac calibration both the polyethylene chambers and the graphite chambers 
were irradiated to a graphite dose level of B rads (0. 01 B J/kg) as measured by the stan- 
dard calibration chamber in a 7. 0-MeV bremsstrahlung beam. This beam contained no 
neutrons because the maximum beam photon energy was below the threshold for y, n re- 
actions. 

Therefore, the readings on the chambers using equations (11) and (12) are 



(13) 


(14) 


where R' and R/, are the scale readings of polyethylene and graphite chambers when 

P L / \ 

irradiated to a dose of B gamma rads (0.01 B J/kg) in graphite, and /E> h J/ D^J is the 

ratio of gamma dose deposited in water to gamma dose deposited in graphite. This con- 
version factor is required since the response coefficients e and K are defined as the 
response per rad (J/kg) gamma dose in water. The response coefficients e and K are 
therefore: 


f 



and 



(15) 


( 16 ) 
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The average value of the polyethylene chamber response coefficient e is computed 
to be e = 1.26 in table V using = 0.9009. The full scale rad (J/kg) val- 

ues for the polyethylene chambers shown in table V were determined in the calibration of 
these chambers against the water-filled calorimeter in TRIGA reactor which is discussed 
later. 

The graphite chambers were absolutely calibrated by this irradiation and the full 
scale reading of each chamber determined by dividing the standard calibrating chamber 
measured dose B by the decimal fraction of full scale read on the chamber . Thus R£ 
was set equal to B in equation (14) and therefore the response coefficient 
K = D c r /D R Q y = 0.9009. 

In the calibration in the TRIGA reactor, the reading of each chamber, R^ and R£ 
was compared against the total radiation deposition C, in rads (J/kg) measured by the 
temperature rise of the water mass in the calorimeter. Thus, 

BJJ = eG" + 5N” 

R ” = KG” + i/N" 

and 

C = G” + N" (19) 

where G” is the rad (J/kg) gamma dose deposited in the water at the position of mea- 
surement in the TRIGA reactor and N" is the rad (J/kg) neutron dose deposited in the 
water at the position of measurement in the TRIGA reactor. 

The polyethylene chambers were absolutely calibrated by this irradiation and the full 
scale reading of each chamber determined by dividing the calorimeter measured dose C 
by the decimal fraction of the full scale reading on the chamber. Thus, R^ was set equal 
to C for this irradiation. The results of this calibration are given in table VIII. 

There are, then, three unknowns in equations (17) to (19), namely, 6, G”, and N”. 
The response of the graphite chamber per rad (J/kg) gamma dose in water v was calcu- 
lated as the ratio of the neutron dose per gram of carbon dioxide to neutron dose per gram 
of water in the section Correction for Neutron Sensitivity of Graphite Chambers. 

By simultaneous solution, the values for these unknowns for the calibration are ob- 
tained as 

R« - vC 

G" = — (20) 

K - v 


(17) 

(18) 
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R£ - KC 


( 21 ) 


N" = 

v - K 


and 


_ R" - eG" C(v - K - Tv) + TRX 

6 = -2 = k (22) 

N" R£ - KC 


The average value calculated for 6 using the values in table V was 6 = 0.861. 
Thus, the general equations with the response coefficients inserted are 


Rp = 1.26G + 0.861N 


(23) 


R c = 0.901G + 0.0634N 


(24) 


These equations were then used to determine the absolute values of G and N from the 
measurements in the TRIGA reactor and in the TWMR critical assembly. 

Methods for obtaining gamma and total dose in TWMR critical assembly . - The meth- 
od for obtaining the absolute gamma dose was explained in the section Correction for Neu- 
tron Sensitivity of Graphite Chambers to require the solution of equation (1) for KG and 
the substitution of the spectrum index S s G/(G + N) in the equation to obtain 


KG = R c 




( 2 ) 


The unknown factor in this equation is the average value of the spectral index S for 
the TWMR critical assembly. This was obtained by solving equations (23) and (24) for N 
and G and then computing the value of S = G/(G + N) at 23 polyethylene chamber core 
locations and then computing the core average. The equations for G, N, and S as ob- 
tained from equations (23) and (24) are 


N = 1 . 29R - 1.81R,-, 
P c 

(25) 

G = 1 . 24R r , - 0. 0911R 
^ P 

(26) 

G 1.24R C - 0.0911Rp 

(27) 

N + G 1.20R^ - 0.570RP 
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TABLE IX. - DATA MEASURED WITH POLYETHYLENE CHAMBERS IN TWMR AND 


CALCULATION OF AVERAGE SPECTRUM INDEX S 


Poison tube or 

Distance from core 

Chamber reading with 

Interpolated 

Spectrum 

fuel element 

bottom plate to 

run normalization 

values of R c 

index, 


center of chamber, 
cm 

included, R p 



G/(G + N) 


rads 

J/kg 

rads 

JAg 


G-2 

67.0 

3000 

30.0 

935 

9.3 

0.289 

G-4 

a 32. 4 



— 

— 

— 



67.0 

4540 

45.4 

1210 

12.1 

0.228 

g-4 

11.4 

2345 

23.4 

810 

8.1 

0.336 


32.4 

3150 

31.5 

930 

9.3 

.267 


67.0 

3070 

30.7 

930 

9.3 

.278 


80.0 

2670 

26.7 

795 

7.9 

.270 


105.3 

1105 

11.1 

340 

3.4 

.284 

g-8 

11.4 

3223 

32.2 

1000 

10.0 

0.268 


32.4 

4380 

43.8 

1355 

13.5 

.286 


80.0 

3360 

33.6 

1000 

10.0 

.270 


105.3 

1490 

14.9 

445 

4.5 

.271 

g-12 

a ll. 4 

3370 

33.7 

1260 

12.6 



a 32. 4 

5030 

50.3 

1765 

17.7 



a 67.0 

— 

— 

— 

— 



a 80. 0 

4040 

40.4 

1230 

12.3 



a 105. 3 

1470 

14.7 

510 

5.1 


h-8 

a ll. 4 



— 

— 

— 



32.4 

4210 

42.1 

1360 

13.6 

0.304 


a 67. 0 

— 

— 

— 

— 



80.0 

3660 

36.6 

1180 

11.8 

.304 


105.3 

1360 

13.6 

480 

4.8 

.347 

j-4 

11.4 

2405 

24.1 

860 

8.6 

0.354 


32.4 

3630 

36.3 

1210 

12.1 

.319 


67.0 

3510 

35.1 

1050 

10.5 

.272 


80.0 

2730 

27.3 

840 

8.4 

.283 


105.3 

1305 

13.1 

350 

3.5 

.231 

k-1 

11.4 

2285 

22.9 

760 

7.6 

0.318 


32.4 

3180 

31.8 

1010 

10.1 

.297 


67.0 

2890 

28.9 

880 

8.8 

.279 


105.3 

1034 

10.3 

345 

3.5 

.320 

g-16 

23.5 

4540 

45.4 

— 

— 



23.5 

4620 

46.2 

— 

— 



23.5 

4530 

45.3 

— 

— 



23.5 

4620 

46.2 

— 

— 


Average 






0. 290±0. 026 


a Suspected defective ionization chamber. 
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Since the computations of S are made at the core locations of the polyethylene 
chambers the values of R can be directly obtained from table IX. The values of R r 

P Ks 

were obtained by interpolation of the plotted measurements of R^ given in table VI. 

Table IX shows the computation of S at the 23 core locations and computation of the 
average value S. The average value of S in the TWMR is S = 0.290, which is then sub- 
stituted into equation (2) and the rad (J/kg) gamma dose in graphite KG computed using 
equation (6) 


KG = 0.853R C 


COMPARISON OF CALCULATION AND MEASUREMENT OF 
ABSOLUTE NEUTRON DOSE 


The neutron dose was computed from the experimental data given in table IX using 
equation (25). The measured results are compared with the calculated results in 
table Vn. 

The calculated values of the neutron dose are consistently low in comparison with the 
measurements which are considered reliable with a deviation from the measurement be- 
tween -11 and -28 percent in 15 out of the 23 locations investigated. 

The error in the neutron dose measurements was estimated on the basis of the un- 
certainty in the polyethylene and graphite dosimeter measurements and on the uncertainty 
of the absolute power measurement. Using equation (25), the uncertainties in the cham- 
ber readings are related to the uncertainty in the neutron dose AN by 


AN = 




(28) 


where the constants are determined using typical values of N = 3200, Rp = 4380, and 
R r = 1355. The uncertainty in R consists of a ±2 percent uncertainty in the reading of 
the polyethylene dosimeter, an ±8 percent uncertainty in the absolute calibration of the 
calorimeter, and a ±3 percent uncertainty in the time at which the TWMR critical assem- 
bly was at power. Combination of these uncertainties results in a ±9 percent error in the 

value of R . 

P 

The comparable error in R^-, is ±6 percent and was discussed in the section Error 
Analysis of Gamma Dose Measurements. 

The error in N obtained from equation (28) is ±17 percent. If the ±7 percent un- 
certainty in the absolute power level is incorporated, the total error in the measured neu- 
tron dose is ±18 percent. Figure 17 shows typical comparisons between the calculated 
and measured neutron dose distribution. 
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COMPARISON OF CALCULATION AND MEASUREMENT OF TOTAL DOSE 


Because of limitations on core storage, ATHENA calculations were made only at the 
positions of the graphite chambers and so a direct total dose calculated at the positions of 
the polyethylene chambers was not obtained. However, since good agreement was found 
between the measured and calculated gamma dose as discussed in the preceding section 
and since the gamma dose contributes only about 30 percent of the total dose, only a small 
error would result if the measured gamma dose were used in place of the calculated dose. 
Thus, the calculated total radiation dose was computed as the sum of measured gamma 
dose and the calculated neutron dose. The calculated total dose computed in this way is 
compared with the measured total dose in table X. The calculated values of the total dose 
are consistently smaller than the measured values. The calculations underestimated the 
measurement by less than 17 percent in 17 out of the 23 locations investigated. This dif- 
ference is caused by the underestimation of the neutron dose discussed earlier. 


MEASUREMENT OF ABSOLUTE POWER LEVEL 

The measurement of the absolute power level was based on the measured critical 
assembly relative power distributions and a calculated conversion factor relating the sub- 
cadmium activity of gold foils to the fission power at the center of the critical assembly. 
The radial power distributions were measured in detail by counting the fuel element fis- 
sion product activity in one symmetrical one twelfth sector at the midplane of the critical 
assembly (see fig. 10). The relative axial power distributions were measured at 24 heights 
at three different fuel element locations also by counting the fission product activity. 

The axial and radial power distributions measured were found to be independent so 
that the local power density to average core power density at the point (r, z) was simply 
the product of the relative radial power density at the radius r times the relative axial 
power density at the height z divided by product of the volume integrated relative radial 
and axial distributions. Table n gives the measured relative radial power density at the 
midplane of the TWMR critical assembly and figure 6 shows the measured axial power 
distribution in the center fuel element. 

The conversion factor relating the subcadmium activity to the power in the center of 
the assembly was evaluated using GAMBLE (ref. 9) computed fluxes and transport- 
theory-determined disadvantage factors appropriate to the position of the foils adjacent 
to center fuel element at the midplane of the core. 

The absolute core power for a reactor run was then determined by dividing the mea- 
sured absolute saturated subcadmium activity by the product of the calculated conversion 
factor and the measured local to average core power at the center of the assembly. 
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TABLE X. - COMPARISON OF CALCULATED AND MEASURED TOTAL 


DOSE IN TWMR CRITICAL ASSEMBLY 


Control tube 
or fuel 
element 

Distance from core 
bottom plate to 
center of chamber, 
cm 

Total heating dose in water 

Percent 

deviation 

Calculated 

Measured 

rads 

J/kg 

rads 

J/kg 

G-2 

67.0 

2950 

29.5 

3060 

30.6 

-4 

G-4 

a 32. 4 

3750 

37.5 

3340 

33.4 



67.0 

3620 

36.2 

4800 

48.0 


g-4 

11.4 

2300 

23.0 

2330 

23.3 

-1 


32.4 

2790 

27.9 

3270 

32.7 

-15 


67.0 

2670 

26.7 

3160 

31.6 

-15 


80.0 

2350 

23.5 

2750 

27.5 

-15 


105.3 

952 

9.5 

1130 

11.3 

-16 

g-8 

11.4 

2880 

28.8 

3300 


mm 


32.4 

3710 

37.1 

4480 




80.0 

3400 

34.0 

3470 


419 


105.3 

1220 

12.2 



-21 

g-12 

a ll. 4 

3350 

33.5 

3250 

32.5 

+3 

* 

a 32.4 

4370 

43.7 

5080 


-14 


a 67. 0 

3950 

39.5 

3050 


+30 


a 80. 0 

3250 

32.5 

4130 

41.3 

-21 


a 105. 3 

1370 

13.7 

1450 

14.5 

-6 

h-8 

a ll. 4 

3040 

30.4 

2180 

21.8 

+39 


32.4 

3750 

37.5 

4250 

42.5 

-11 


a 67. 0 

3550 

35.5 

2740 

27.4 

+30 


80.0 

3020 

30.2 

3700 


-18 


105.3 

1260 

12.6 

1310 

13.1 ’ 

-3 

j-4 

11.4 

2430 

24.3 

2360 

23.6 

+3 


32.4 

3170 

31.7 

3840 

36.4 

-13 


67.0 

2880 

28.8 

3620 

36.2 

-20 


80.0 

2360 

23.6 

2800 

28.0 

-16 


105.3 

994 

9.94 

1380 

13.8 

-28 

k-1 

11.4 

2300 

23.0 

2290 

22.9 

+0 


32.4 

2920 

29.2 

3220 

32.2 

-9 


67.0 

2670 

26.7 

2960 

29.8 

-7 


105.3 

974 

9.74 

1050 

10.5 

-7 


a Suspected defective ionization chamber. 









An absolute power level of 84.2 watts was measured. An error of ±7 percent was 
assigned to this measurement based on the following estimates of individual errors: 

±2 percent in the measurement of the subcadmium ratio of gold, ±0. 2 percent in the acti- 
vation cross section of gold, ±4 percent in the calculated conversion factor, and ±5 per- 
cent in the measured local to average power distribution. 


CONCLUSIONS 

The mixed gamma and neutron radiation dose can be measured accurately in a reac- 
tor using graphite wall, carbon dioxide filled, and polyethylene wall ethylene filled ioniza- 
tion chambers. The mixed radiation dose can be partitioned to yield separately the value 
of the gamma and the neutron dose if the ionization chambers are absolutely calibrated 
both in a pure gamma flux and in the mixed spectrum of a reactor. In the work reported 
herein, the ionization chambers were absolutely calibrated in a bremsstrahlung spectrum 
and against a water filled calorimeter in the reactor. In addition, the modest neutron 
sensitivity of the graphite chambers must be determined. In this work, the neutron sen- 
sitivity was calculated. 

The agreement between the ATHENA calculated and the measured gamma dose was 
excellent with a deviation within 10 percent in 10 out of the 16 locations compared. The 
absolute calculated neutron dose was uniformly low with a deviation from the measure- 
ment between -11 percent to -28 percent in 15 out of the 23 locations compared. The 
cause of this underestimation is not known but may be due to the refined first collision 
approach used in the calculation. The calculated combined total gamma and neutron dose 
underestimated the measurement by less than 17 percent in 17 out of the 23 locations in- 
vestigated. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 22, 1967, 

120-27-06-18-22. 
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APPENDIX - SYMBOLS 


A(r ) 

B 

C 

D 

d h ? o/ d c 

E 

AE 

E i 

E 1 

F A< Z > 

f r (D 

f 

G 

K 


M 

N 

AN 

*1 

n 

R C 

R P 


r 


normalizing factor relating average dose at r to average core dose 
irradiation dose level of bremsstrahlung beam 
heat deposition in water mass of calorimeter rad; J/kg 
dose, rad; J/kg 

ratio of gamma dose deposited in water to gamma dose in graphite 

total dose deposited in 1 gram of material per unit time from fast neutrons 
at position r, rad; J/kg 

energy, eV 

energy interval, eV 

energy at midpoint of each interval, eV 
lower energy limit of fast neutron range 

axial factor corresponding to power at distance z above core bottom plate 

radial factor corresponding to relative power density at radius R 

average fraction of energy imparted to recoiling atom per elastic collision 

o 

with neutron equal to 2M/ (M + 1) 

gamma dose deposited in water at position of measurement, rad; J/kg 

ratio of gamma dose in graphite to gamma dose in water; also, gamma dose 
coefficient or dose reading of graphite chamber per rad (J/kg) gamma 
dose in water 

atomic mass 

fast neutron dose deposited to water at position of measurement, rad; J/kg 
uncertainty in neutron dose 

j. i_ 

number density of i in isotope, atoms/g 
neutron 

graphite chamber reading for mixed radiation dose of G rad gamma dose 
plus N rad neutron dose in water, rad; J/kg 

polyethylene chamber reading for mixed radiation dose of G rad; J/kg 
gamma dose plus N rad neutron dose measured in water, rad; J/kg 

radius 
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r position vector 

S spectrum index G/ (G + N) 

u lethargy 

Au lethargy width equal to 0.1 from 14.9 MeV to 86.5 keV and 0.25 from 86.5keV 

to E l 

z coordinate 

5 neutron dose coefficient or dose reading of polyethylene chamber per rad 

(J/kg) neutron dose in water 

e gamma response coefficient or dose reading of polyethylene chamber per rad 

(J/kg) gamma dose in water 

v ratio of fast neutron dose (rads or JAg) in carbon dioxide to fast neutron dose 

(rads or JAg) in water; dose reading of graphite chamber per rad gamma 
dose in water or neutron response coefficient 

2 

ct_ elastic scattering cross section, b; m 

s 

cp(E) flux per unit energy interval about energy E 

( p(E , r ) neutron flux per unit energy about energy E at position r 

(p (U) neutron flux per unit lethargy 

Subscripts: 

C carbon or graphite 

CC>2 carbon dioxide 

H hydrogen 

H2O water 

i fine group 

j fine group 

o oxygen 

’ in LINAC calibration 

” in TRIGA reactor 
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" The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space . The Administration 
* shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

* — National Aeronautics and Space Act of 1958 
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